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We present non-perturbative renormalization constants of fermionic bilinears on the lattice in the 
quenched approximation at (3 = 6.1 using an overlap fermion action with hypercubic(HYP)- 
blocked links. We consider the effects of the exact zero modes of the Dirac operator and find they 
are important in calculating the renormalization constants of the scalar and pseudoscalar density. 
The results are given in the RT and MS schemes and compared to the perturbative calculations. 
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This paper describes the computation of matching factors for converting lattice calculations of matrix elements of 
currents to the corresponding values measured in a continuum MS scheme. The lattice action uses overlap fermions pj 
and HYP-blocked links 0. 

^v^j , We use the nonperturbative methodology introduced in Ref. |3( . The proposed renormalization scheme is one which 
can be implemented not only in lattice Monte Carlo simulation but also in continuum perturbation theory. Thus, 
the conversion of lattice results to a more conventional scheme such as MS is possible. In this scheme, the matrix 
clement of a bilinear quark operator Or = TpTij) between quark fields at certain momentum p 2 = /i 2 is computed and 
matched to the corresponding tree level matrix element, i.e. the renormalization condition is 

Zr< P |Pr|p)| P »= M »=<p|Or|p) tree . (1) 
. Here L can be any combination of Dirac matrices. This method is supposed to work when /i satisfies 



o 
o 
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Aqcd < A* < 1/a- (2) 



The discretization effects are under control if the renormalization scale [i is much smaller than the lattice cut off 1/a. 
'. h-QCD <C \i guarantees that the non-perturbative effects are ignorable. 
• There have been many calculations using this scheme. Ref. Q used improved Wilson fermions, Ref. £| used both 
the Wilson and the tree level improved SW-Clover fermion action in the quenched approximation, Ref. [f| worked with 
standard Wilson fermions(r = 1) in the quenched approximation, Ref. [6j used chirally improved lattice fermions in 
the quenched approximation. Here, we use overlap fermions Q in the quenched approximation. Specifically, we work 
Oh, with overlap fermions built from a "kernel action" with nearest and next-nearest neighbor fermionic interactions Q 
and hypercubic(HYP)-blocked links 0j. Overlap fermions respect chiral symmetry on the lattice via the Ginsparg- 
Wilson relation Q, while for Wilson-type fermions, the Wilson term breaks the chiral symmetry explicitly. Chirally 
^ ■ improved fermions only obey the Ginsparg- Wilson relation approximately. 

Lattice perturbation theory is probably the most often used way to calculate the renormalization factors. However, 
the convergence of the perturbative series is often not satisfying. To improve the convergence of the series, Lepage 
and Mackenzie proposed a tadpole improved perturbation theory [(| . Nevertheless, lattice perturbation series rarely 
extend beyond the one-loop level, which is an important source of uncertainty in the extraction of physical results. 

One-loop perturbative calculations of the matching coefficients between matrix elements measured in lattice simu- 
lations and their equivalent MS values for the same overlap fermions and HYP-blocked links that we use here were 
done in Ref. |10( . Those perturbative results turned out to be quite close to unity, and they were used in computing 
the Kaon B parameter This work will give a non-perturbative check of the matching coefficients. 

Perturbative results of the matching coefficients for other actions also using HYP-blocked links or similar gauge 
connections were presented in [l2 |. They show the same behavior that the matching coefficients are quite close 
to unity. Our results may be useful to others doing simulations with HYP links, to give an idea how trustworthy 
perturbation theory is. 

In Eq. r can be any combination of Dirac matrices. We will consider the cases T = I, 75, 7^ and 7^75, which we 
will denote S, P, V and A respectively. Chiral symmetry implies several relations between renormalization constants 
for overlap fermions, in particular Zg — Zp and Zy = Z&. 

The paper is organized as follows: In Section[Hj we briefly discuss the non-perturbative renormalization method 0], 
the overlap action and how we deal with the zero modes of the Dirac operator. Numerical results are given in Section 
IIIII The formulas for conversion to A/5 scheme are recapitulated in Section IIVI We will compare our results with 
perturbative calculations in Section and conclude in Section IVII 
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II. METHODOLOGY 

The following is a brief summary of the method from Ref . 0, , a short description of the overlap action we used 
(for a detail description, see Ref. Q) and how we deal with the zero modes of the Dirac operator. For convenience, 
the lattice spacing a is set to be one. 

From Eq. (JTJ, we have 

Zr-^K [{p\0 T \pMO T \p)^ ee } \ p2 ^ 2 = 1. (3) 

Here jk comes from the fact that the trace is over color and spin space. Since 

(p\O r \p) = Z q A r (p), (4) 

we obtain 

12 



Zr = 



^TrfArbXplOrk^-L 



(5) 



1 /2 

Here Z q is the quark field renormalization constant (The bare field = Z q ip) and Ar(p) is the amputated Green 
function 

A r (p)=,S- 1 (p)G r (p)S-\p), (6) 

where S(p) is the quark propagator. Eq. © is the formula we will use to calculate Zr- 

Z q is obtained by comparing the quark propagator to the free lattice propagator (the RT scheme): 



Zf = ±Tr[S(p)D?(p)] 



(7) 



where D°f V {p) is the free lattice overlap Dirac operator. (Our Z q is the inverse of the quark field renormalization 
constant in Ref. 0.) 

The Green function Gr (p) is determined in the following way. 

x,y 

1 N 

x,y i=l 

= jf E (E s,(x\o)e-^ r (^2 SMy)e ip v ^ , (8) 

where N is the number of gauge configurations. Using Si(x\y) = r YnSi(y\x) 1i ~f 5 and Si(p\0) = J2 X Si(x\0)e~ lp x , we 
have 



1 N 

Gr W = n E ^(P|0)r 7 5^(p|0)75. (9) 

i=l 

The quark propagator in momentum space is given by 

1 N 

Si(x\0) is computed on the lattice with a point source 

J2D ov (z,x)Si(x\0) = 5 z , . (11) 
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(In Ref. [f| and Ref. momentum sources were used.) At tree level, (p\Or\p)tree = T. Therefore, every quantity on 
the right hand side of Eq. (JjjJ is known and then we can obtain Z?. For T = 7^,7^75, the index 11 is averaged under 
the trace in Eq. 

We fix the gauge to Landau gauge. Uncertainty due to Gribov copies is not investigated here. It has been discussed 
in Ref. |rj lilj. TbjTll5l|. The effect was found to be negligible in current lattice simulations. 

The overlap action that we use is described with detail in Ref. Q, which uses a "kernel" action with nearest and 
next-nearest neighbor couplings. The massless overlap Dirac operator is 



D{0)= X0 [l + 7T J, (12) 

where z = d(—xo)/xo = (d — xq)/xq and d(m) = d + m is the massive Dirac operator for mass m. The overall 
multiplicative factor of xq is a useful convention so that when D(0) is expanded for small d, D ss d. 
The massive overlap Dirac operator is defined as 

(vn \ 
1 )D(0) + m. (13) 
2x J 

In a background gauge field carrying a topological charge Q, D(0) will have \Q\ pairs of real eigenmodes with 
eigenvalues and 2xo- In computing propagators, it is convenient to clip out the eigenmode with real eigenvalue 2xo, 
and to define the subtracted propagator as 



Dim)- 1 = — 



Dim)- 1 - — 
2x 



(14) 



2x 

This also converts local currents into order a 2 improved operators ^(|. Then the free lattice overlap Dirac operator 
D^ip) used in Eq. (J7J) is just £>(m) in the momentum space. 

The HYP-blocked links are constructed in three steps 0. The parameters ot\, oti and in our simulation have 
the favored values of Ref. Q: 0.75, 0.6 and 0.3 respectively. 

A finite volume artifact we encounter in this quenched simulation is the presence of exact zero modes of the Dirac 
operator. The zero mode contribution (with positive chirality) in the propagator Si(p\0) on a configuration with 
Q 7^ takes the form 

L(\MpMMp)\ ;U So 

in 75 diagonal basis. Here \<po(j>)) is the Fourier transform of the zero mode wave function \<j)o{x)). Since the zero 
modes are localized in space, |<^o(p)) will peak at low p. These zero modes do not resemble free field modes. Implicit 
in the RT scheme analysis is the idea that at big lattice propagators resemble continuum ones. Zero modes clearly 
do not. In Section ITTT1 we will find zero modes make a large contribution to Zg and Zp. We believe this is because 
our lattice is not large. 

The following little parametrization illustrates our expectations of the effects of zero modes: Si(p\Q) is the sum of 
the zero mode contribution and the non-zero mode contribution S n , 

Siip\0) = ^S + S n . (16) 

Therefore in Eq. © 

^10^75^10)75 = ^S T l5 Sl l5 + ^(5^75^75 + S n T l5 S^ 5 ) + S n T~/ 5 Sty 5 , (17) 
and then Gr{p) can be written in the form 

G T ip) = X [ G 2 + -G 1 + G , (18) 

where the subscript counts the number of zero modes: Go contains no zero mode contribution. The quark propagator 
averaged over all configurations and its inverse, if expanded for small m, are 

S(p) = — s + s n , 

m 

S^ip) = mS^ 1 - m'S^SnS^ 1 + ■ • ■ • (19) 
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Thus the amputated Green function 



A r (p) = S-^GrWS- 1 ® 
1, ,VVl 



— So + Sri — 5-G*2 H Gi + Go — So + <5 ra 

m J \m m 

Sq 1 G2Sq 1 when m is small. 



1 

m 



(20) 



So, if the zero modes affect our calculation of Z-p (Eq. (JSJ), the effect should be evident at small momentum and small 
quark mass. Unfortunately, for us, fi = 2 GeV, where we will match, is rather small momentum. 

We examine two solutions to those zero modes. One solution is to explicitly subtract the contribution of the zero 
modes in the quark propagator. We would prefer not to do this. The other solution is to use combination of scalar 
and pseudoscalar densities or combination of vector and axial vector currents so that the zero mode contributions are 
suppressed in Zr- We can do the latter because the overlap fermion respects chiral symmetry on the lattice. 

For scalar and pseudoscalar, we can use Eq. to rewrite Eq. as 



and 



If Zsp = Zs = Zp, then we have 



and thus 



Eq. I© and Eq. I© give us 



Z S Z q Al(p 2 = ^)=I 

Z P Z q Ay 5 (p 2 = n 2 ) = 75. 
Z SP Z q (A I ±A j5 ) = I± l5 , 



Z SP Z g -T I (A I ±A 75 ) = l. 



A/±A 75 =S-\ P ) 



1 N 

-£>(p|0)(/± 75)755/(^10)75 
i=i 



(21) 
(22) 
(23) 
(24) 

(25) 



For zero modes with positive chirality, Si(p\0)(I— 75) = 0, while for zero modes with negative chirality, £j(p|0)(i+75) = 
0. Therefore, the zero mode contribution in the Green function combination Gs{p) ± Gp(p) are removed. We can 
use Eq. I|24|l to obtain a Zsp which has a suppressed zero mode contribution (Note that S~ 1 (p) still contains zero 
modes) . 

Similarly, for vector and axial vector currents, Eq. gives 



and 



Letting Zva = Zy = Za, we obtain 



ZyZqA^^p 2 = y 2 ) = 7 M , 
ZAZ q A 1 ^(P 2 = M 2 ) = 7m75- 



j^Z VA Z q (A 7fi ± A 7ji75 ) = I ± 75 
and subsequently (after the index /i is averaged) 



Z VA Z q —Tr [(A 7fi ± A 7ti75 ) 7M ] = 1. 



Similar to Eq. H25(l . we have 



1 N 

-^5(p|0) 7 , 1 (/±75)755 t (p|0) 7 5 



s-'ip). 



(26) 
(27) 
(28) 

(29) 
(30) 



For zero modes with positive chirality, Si(p\0)-f tI (I + 75) = Si(p\0)(I — 75)7^ = 0. For those with negative chirality, 
5(^10)7^(7 — 75) = Si(p\0)(I + 75)7^ = 0. Therefore, Z VA calculated from Eq. has a suppressed zero mode 
contribution. 



5 



1,0 



0.9 



0.8 



0.7 



0.0 0.5 



1,0 



am q =0.020 



0.9 



O full propagator 

n zero modes subtracted 



0.8 



1.0 1.5 
afj, 



0.7 

2.0 2.5 0.0 0.5 



am q =0.070 



O full propagator 

n zero modes subtracted 



1.0 1.5 
a/j. 



2.0 2.5 



FIG. 1: Zq J vs. afi for bare quark mass am, = 0.020 and 0.070. Z^ 1 obtained from the full propagator (diamond) is 
compared with that obtained from the propagator with zero mode contribution subtracted (square). 



III. NUMERICAL RESULTS 



The data set that we use contains 40 gauge configurations in the quenched approximation with the Wilson gauge 
action. The lattice size is 16 4 and the gauge coupling /3 = 6.1. The bare quark masses in lattice units are am q = 0.015, 
0.020, 0.025, 0.035, 0.050 and 0.070. The lattice spacing a is (a): 0.08 fm determined from the interpolation formula 
of Ref. 01 using the Sommcr parameter or (b): 0.09 fm from the measured rho mass. Therefore, a/i = 0.811 or 
a\x = 0.913 corresponds to fi — 2 GeV accordingly. In the following analysis, the statistical errors are obtained by a 
Jackknifc average with one configuration removed each time. 



Z* 1 ' and Z* 1 ' 



The quark field renormalization constant Z^ 1 is calculated with Eq. Q . The results for two examples of bare quark 
masses are shown in Fig. ^ The comparison between Z^ 1 obtained from the full propagator and the propagator with 
zero mode subtracted is also shown in the same graph. There is no difference within error bars. 

The full lattice quark propagator takes the form 

s ^ = ~- — r¥rw( y ( 31 ) 
h ■q{p) + M{p) 

Here q(p) is the kinematic momentum depending on the lattice quark action one uses. At large momentum p, because 
of asymptotic freedom the propagator should go back to the free quark propagator, i.e. Z(p) — ► 1 and M(p) goes 
to the bare quark mass. Fig. |2Jshows (l/12) r Tr(5 _1 (p)) versus ap for two examples of bare quark masses with S(p) 
determined from Eq. I|10|) . Results from the full propagator and from the propagator with zero mode contribution 
subtracted are compared in the graph. As is expected, (l/12)Tr(S~ 1 {p)) approaches the bare quark mass at large 
momentum. Apparently, only at small momentum and small quark mass does the zero mode contribution make a 
difference. 

If we define a renormalizcd quark mass m(fi) by 

m(fjb) = Z m (n)mo, (32) 
where too is the bare quark mass, then Z m {\i) is fixed in the RF scheme by 



p 2 =ii 2 



(33) 
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FIG. 2: (l/12)Tr(S _1 (p)) vs. ap for bare quark masses am q — 0.020 and 0.070. Results from the full propagator (diamond) 
and from the propagator with zero mode contribution subtracted (square) are compared. Zero mode contribution is important 
only at small momentum and small quark mass. 
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FIG. 3: (Z^l ) 1 for bare quark masses am q = 0.020 and 0.070. Diamonds are from the full propagator and squares from the 
propagator with zero mode contribution subtracted. 

At finite quark masses, the renormalization conditions of RT scheme are compatible with the Ward identities 0, 0] 
at large /it 2 , therefore we expect = Zg 1 at large pi 2 . The numerical results of (Z^ 1 ) _1 are shown in Fig. The 
error bar at small quark mass is large. A zero mode contribution is visible only at low momentum. We will compare 
(Z™')- 1 with Z s later. 



B. Z§'' and ZjF' 



The results for Zs and Zp are shown in Fig. At low quark mass and momentum region, the zero mode subtracted 
propagators give different values of Zs and Zp (we will label them as Zg Z and Zp z in the following) from those 
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FIG. 4: Zs and Zp from the full propagators and the zero modes subtracted (labeled as Zg Z and Zp Z ) propagators. The left 
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FIG. 5: Zsp obtained by using Eq. 12 It comparing with Zs and Zp. The left is for quark mass am q — 0.020, the right 
am q = 0.070. Q is the topological charge. 



obtained with the full propagators. The pseudoscalar density couples to the Goldstone boson channel but the coupling 
is suppressed at large /i j^]. Therefore we see a difference between Zs and Zp at small /i, but no difference at large 

As was discussed in Section [n] we can also use Eq. I|24f) to suppress the zero modes and obtain Zgp. Fig. [S] shows 
the results of Zgp comparing with Z$ and Zp. The Zsp's from configurations with topological charge Q > and 
Q < are quite close to each other. At small quark mass and small fj,, Z$p is apparently different from Z$ and Zp 
as can be seen in the graph for am q = 0.020 in Fig. EI Zsp agrees with Zg Z and Zp Z obtained from the zero mode 
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FIG. 6: (Z* 1 ')' 1 compared with Z^ 1 ' . At large /i, (Z* 1 ')' 1 
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Zg 1 is well satisfied. 
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FIG. 7: Linear extrapolation of Zs to the chiral limit. The left is for Zg obtained from full propagators. The data do not 
support a linear extrapolation. We did it anyway just for comparison. The right is for Zs obtained from zero mode subtracted 
propagators (labeled as Zg Z ). At small quark mass, Zs and Zg Z are very different. 



subtracted propagators in Fig. fj] for am q = 0.020. For large quark mass, for example am q — 0.070, and small fi, Z$p 
is close to Zs but very different from Zp. Z$p from Eq. (|24f> is still contaminated by the coupling to the Goldstone 
boson. Thus, this means that suppressing the zero modes can suppress the coupling to the Goldstone boson. We 
do not see this behavior in Fig. [J] for am q = 0.070. The Zp z from the zero mode subtracted propagator is very 
close to the Zp from the full propagator. However, we should notice that subtracting zero modes directly from the 
propagators amounts to a modification of the quenched theory. To further investigate the zero modes, a quenched 
artifact, simulations with dynamical fcrmions are necessary. 

The comparison of (Z^ 1 )~ 1 with Zg 1 is given in Fig.El We see a good agreement between (Z^ 1 )~ 1 and Zg 1 at 
large [i as expected. 
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FIG. 8: Extrapolation of Zp to the chiral limit using Eq. 1341 . Zp in the left graph is obtained by using the full propagator, 
while Zp Z in the right graph is obtained by using the propagator with zero modes subtracted. Both are fits after the last 
configuration is dropped during the Jackknife average process. 



In Fig. d O and El Zg, Z^ z , Zp, Zp Z and Zsp(average from configurations with Q > and Q < 0) are plotted 
versus the quark mass at afj, — 0.708, 1.057 and 1.583 along with the extrapolation to the chiral limit. A linear 
fit is used for Z§ and Zg Z . The data for Zs which is obtained from the full propagator do not support a linear 
extrapolation as is shown in Fig. For Zp, we see similar behavior as was seen in 0, 0, [53, EJ since the 
pseudoscalar density couples to the Goldstone boson channel. As in |2lJ , we use 



] +B( j i 2 ) + C( f i 2 )(am q ) (34) 



Zp(fi 2 ,m) am q 

to fit Zp and then remove the pole term A(/j, 2 ) / am q to obtain Zp P = B(fi 2 )~ 1 in the chiral limit. The fit is good as 
can be seen in Fig. |SJ Fig. [S] shows one example of the fit when we drop one configuration and obtain Zp with the 
rest configurations during the Jackknife average process. 

We also use Eq. <|34ll to extrapolate Z$p to the chiral limit to obtain Zg P = B(n 2 )~ l . The fit is shown in Fig. 
Values of Z s , Z^ p , Zg z , (Z P JZ ) NP and Zg p in the RT scheme in the chiral limit arc listed in Table [Q In the 
table, the superscripts (a) and (6) indicate the two different ways of determining the lattice spacing. We could not get 
Zs for the smallest momentum because the signal is noisy in our data. Zs and Zp p arc different from each other at 
small fi so that the MS values of them differ from each other. The MS values are obtained by using the conversion 
formulas in Section Hvl and a linear interpolation from the two closest \x value of the data. After subtracting the zero 
modes from the quark propagators, we get Z$ z and (Zp Z ) NP . They are in good agreement with each other, which 
is expected from chiral symmetry, but are very different from Zs and Zp P . The other way of suppressing the zero 
modes gives us Z$p. Since the zero mode effects arc much more apparent in Zs and Zp, we will use Zg Z , (Zp Z ) NP 
and Zg P to compare with perturbative calculations. 



C. Z^ 1 ' and Zf' 

In Fig. 1101 the renormalization constants of the vector current Zy and axial vector current Za are given. For clarity, 
at each value of the momentum afi, the x-positions of Zy and Za are shifted a little in the graph. As is shown in 
the graph, Zy and Za are independent of the scale at large /i (a/j, > 0.7), and Zy = Za within statistical errors. At 
low ii, Za is bigger than Zy . We think it is because the axial vector current is coupled to the Goldstone boson. Zy Z 
and Zj[ z obtained from quark propagators with zero modes subtracted are compared with Zy and Za in Fig. 1111 
Apparently, zero modes have little effect on Zy and Za- Zy a obtained from Eq. (|29|) is shown in Fig. It agrees 
with the average of Zy and Za- This confirms that zero mode contribution doesn't matter in the computation of Zy 
and Z a- 
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FIG. 9: Extrapolation of Zsp to the chiral limit using Eq. 1341 . Zsp is the average from configurations with Q > and Q < 0. 



0.439 
0.708 
0.982 
1.057 
1.194 
1.583 
2.050 


M (a) (GeV) 
1.08 
1.75 
2.42 
2.61 
2.95 
3.91 
5.06 


/i (b) (GeV) 
0.96 
1.55 
2.15 
2.32 
2.62 
3.47 
4.49 


Zs 

0.85(4) 
0.89(3) 
0.86(2) 
0.89(2) 
0.90(1) 
0.896(7) 


r,NP 

Zj P 

1.31(6) 
0.97(1) 
0.94(2) 
0.887(9) 
0.92(1) 
0.906(4) 
0.901(2) 


0.43(2) 
0.64(1) 
0.753(9) 
0.764(8) 
0.805(5) 
0.860(5) 
0.868(6) 


{Z N Z) NP 

0.455(2) 
0.631(1) 
0.755(1) 
0.7685(9) 
0.817(1) 
0.8611(3) 
0.8725(2) 


r,NP 
^SP 

0.642(2) 
0.732(2) 
0.828(2) 
0.7972(9) 
0.8489(9) 
0.8874(6) 
0.8826(4) 


MS (a> 
MS {b) 


2 GeV 


2 GeV 


1.01(2) 
1.01(1) 


1.12(1) 
1.091(9) 


0.79(5) 
0.83(3) 


0.79(5) 
0.83(4) 


0.89(4) 
0.93(3) 



TABLE I: Values of Z s , Z^ p , Zg z , (Z^ Z ) NP and Zg£ in the RT scheme in the chiral limit. The MS(2 GeV) value is 
obtained from a linear interpolation from the two closest values of the data. The lattice spacing is (a): 0.08 fm from Sommer 
parameter or (fe): 0.09 fm from the measured rho mass. Correspondingly, we get two MS (2 GeV) values. Zs and Zp P contains 
zero modes and seen to be quite different from the Z's which contains no zero modes. 

The RT scheme values of Zy 1 Za : Zy Z and Z^ z in the chiral limit are given in Table ITT1 The superscripts (a) 
and (6) in the table indicate the two different ways of determining the lattice spacing. The MS values are obtained 
by using the conversion formulas in Section IIVI Zy = Za is very well satisfied as expected since the overlap fermion 
respects chiral symmetry on the lattice. The linear extrapolation to the chiral limit is shown in Fig. 1131 

IV. CONVERSION TO MS 

The ratio Zy is (/s 2 )/Zp l (/i 2 ) connects the MS scheme to the RF scheme is computed by continuum perturbation 

theory. There is a need to determine the coupling constant a s (^t) in the ratio. We obtain a^ s (fi) by first measuring 
the trace of the plaquette operator U p i aq (the lxl Wilson loop), which can give us a^(3.41/a) Then aAy and 

aAjjg are calculated. Finally a^ s {y) is determined by 

(aF^(M)r 1 - A) IhWAms) 2 + (A/A)) lnlnOi/A^) 2 , (35) 

where (3q = 11 /An and fi\ = 102/167T 2 for the quenched approximation. If the lattice spacing a = 0.08 fm from the 
Sommer parameter, we find af^Ou = 2 GeV) = af s (0.811/a) = 0.2038. If a = 0.09 fm fr om the measured rho mass, 
then a^ip = 2 GeV) = 0^(0.913/0) = 0.1940. 
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1.3 



1.4 



1.2 



1.1 — 



am q = 0.020 



O Zv 



x 0.5#(Z V +Z A ) 



1.3 — 



1.2 



am q =0.070 



1.1 



O Z v 

x 0.5#(Z V +Z A ) 

+ z a 



1.0 — 



1.0 — 
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S./J, B./J, 

FIG. 10: Zv,Za and their average versus a/j, for quark masses am q = 0.020 and 0.070. For clarity, at each value of the 
momentum afj,, the x-positions of Zv and Za are shifted a little. 



1.3 



1.2 — 



i.l 



1.0 



IT 



am q =0.020 



+ z a « Za 



1.4 



1.3 



1.2 



am q =0.070 



o Z v n Z v 

♦ Z A « Z A 



1.1 — 



1.0 



q g I i i i i I i i i i I i i i i I i i i i I i i i i I q g I i i i i I i i i i I i i i i I i i i i I i i i i 
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a/j, a/j, 

FIG. 11: Zy z and Z% z compared with Z v and Za for quark masses am„ = 0.020 and 0.070. 



For the scalar and pseudoscalar, in the Landau gauge and 3-loop order, the conversion ratio is |18|, |2 



Z¥ s Zl IS 16 a s /4291 152C; 



yRI> yRI 1 



= 1 



3 ^ (®s_ 
47T 



3 4tt V 18 3 
/ 3890527 224993C 3 2960C 5 \ f<M 3 
1 324 54 4 



5 )(S) 3+ °^- 



(36) 



where £„ is the Riemann zeta function evaluated at n. Substituting a^ IS (fi = 2 GeV) = 0.2038 or 0.1940 into the 
above equation, wc get Zf~ & '/Zg r = Zf^/Z^ 1 ' = 1 + 0.08650 + 0.04668 + 0.03131 = 1.1645 or 1 + 0.08234 + 0.04230 + 
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1.3 
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1.1 



1.0 



0.9 
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o Z VJ (20 configs with QSO) 
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FIG. 12: Zva and 0.5(Zv + Za) at the smallest quark mass am q — 0.015. 



0.439 
0.708 
0.982 
1.057 
1.194 
1.583 
2.050 


/i (a) (GeV) 
1.08 
1.75 
2.42 
2.61 
2.95 
3.91 
5.06 


0.96 
1.55 
2.15 
2.32 
2.62 
3.47 
4.49 


Z v 
1.11(3) 
1.03(1) 
1.024(9) 
0.998(9) 
1.026(6) 
1.023(5) 
0.992(5) 


Z A 
1.17(4) 
1.03(1) 
1.026(10) 
1.000(9) 
1.029(6) 
1.023(6) 
0.992(5) 


rvNZ 

1.15(2) 
1.04(1) 
1.023(9) 
0.999(8) 
1.029(6) 
1.023(6) 
0.992(5) 


rjNZ 

A A 

1.16(3) 
1.04(1) 
1.022(9) 
0.999(8) 
1.029(6) 
1.023(6) 
0.992(5) 


MS (a} 
MS (b) 


2 GeV 


2 GeV 


1.022(2) 
1.021(1) 


1.023(1) 
1.022(1) 


1.028(6) 
1.022(4) 


1.028(7) 
1.022(4) 



TABLE II: Values of Z v , Z Al Zy and Z% in the RT scheme in the chiral limit. The MS (2 GeV) value is obtained from a 
linear interpolation from the two closest /i values of the data. The lattice spacing is (a): 0.08 fm from Sommer parameter or 
(&): 0.09 fm from the measured rho mass. Correspondingly, we get two MS(2 GeV) values. 



0.02701 = 1.1516. 

For the vector and axial vector, since Zy 1 



'MS 



and the difference between RI scheme and RF scheme is only 



the different definition of the quark field renormalization constants, we have 



zf 



yMS 
Zy 



7 RI 

7 RI' 
Zy 



Z, 



m 1 



IK 



z** 



Zf 1 ' ' jZf s and Zf 1 /Z^ IS were calculated in Rcf. 



zr_ 

Zf 



Zy 1S 



67/0^x2 / 52321 
~6 \4^J ~ \ 72 



Z^/ZMS 

to 3-loop, so we find 
607C 3 



Utt 



(37) 



0(a 4 s 



(38) 



The numerical value at \ i = 2 GeV is (1 + - 0.00294 - 0.00232) = 0.9947 or (1 + - 0.00266 - 0.00200) = 0.9953. 

In Table HI and ITT1 the MS values at /j, — 2 GeV are obtained from linear interpolations between the two closest fi 
values of the data. 
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FIG. 13: Extrapolation of Zy and Za to the chiral limit. 





0^(1.96/0) 


0^(1.52/0) 


af s (1.96/o) af s (1.52/a) 


Z,s,p 


1.010 


1.011 


1.008 1.009 




0^(1.26/0) 


aY (1.46/a) 


af s (1.26/a) af s (1.46/a) 


Zy,A 


0.989 


0.990 


0.991 0.992 



TABLE III: Values of Zs,p and Zv,a at a/Li = 1 for HYP-planar overlap action from perturbative calculation in Ref. [TfT ]. 



V. COMPARISON WITH PERTURBATIVE CALCULATIONS 

The perturbative calculation in Ref. 0] gives the lattice to MS matching factor Zi = 1 + z,a! s (<i»*)/37r at a/i = 1. 
Here i=S, P, V and A for fcrmion bilinears. The values of Zi and the scale q* are given in Table V in Ref. [lOj. 

We may use aY (<?*) run from a^(3.41/a), as is determined from the plaquette Q or a^ /s (<7*) from Eq. (|35|) . The 
results of Zj's are listed in Table ITTT1 The ambiguity in the choice of a s and q* in perturbation theory is small. 
We have to run the result of Z$,p to /i = 2 GeV to compare with our MS(2 GeV) value. We use the two loop 
formula for the running quark mass given in Ref. [2j| (Eq. (4.81)). If the inverse lattice spacing is 1/a = 2.47 GeV 
(a = 0.08 fm) from the Sommer parameter, then we find Z S , P (2 GeV) = 0.975 from Z S , P (2.47 GeV) = 1.009. If the 
inverse lattice spacing is 1/a = 2.19 GeV (a = 0.09 fm) from the measured rho mass, then Zg.p(2 GeV) = 0.995. 
In any case, the value of Z$,p from perturbative calculation is quite close to 1, while our non-perturbative results 
0.79(5)/0.89(4) or 0.83(4)/0.93(3) (see Table [J are not. Thus, perturbative calculation of the matching factors for 
scalar and pseudoscalar density for HYP-planar overlap action seems unreliable. 

Unlike Zs or Zp, Zy and Za are scale independent. We can compare the values of Zy,A in Table ITTll directly with 
our non-perturbative MS(2 GeV) results in Table UTI All are quite close to one (the shift from one is less than 0.03). 
This indicates that we can believe in the perturbative calculations of Zy and Za for the HYP-planar overlap action. 



VI. SUMMARY AND CONCLUSION 



We calculated the renormalization constants of bilinear quark operators non-perturbatively using the HYP-planar 
overlap action with exact chiral symmetry By comparing the results with those from perturbative computation, we 
find that a perturbative calculation is reliable with Zy and Za, but not with Z$ and Zp. The exact zero modes of 
the Dirac operator turn out to be important in calculating Zs and Zp, while not relevant in calculating Zy and Za- 
After subtracting the zero modes from the quark propagator, Zg = Zp is well satisfied. Zy and Za are also in good 
agreement with each other as is expected from the chiral symmetry of the action. We expect that zero modes will be 
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much less important in simulations done with dynamical overlap auarksp^. 

The perturbative result that actions using HYP-blockcd links have matching factors quite close to unity is confirmed 
for vector and axial vector currents with our HYP-planar overlap action. This does not appear to be the case for the 
scalar and pseudoscalar densities. 
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